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Multifractal analysis of multiplicative random cascades is revis- 
ed , ited within the framework of mixed asymptotics. In this new frame- 
work, statistics are estimated over a sample which size increases as 
the resolution scale (or the sampling period) becomes finer. This al- 
^SJ ■ lows one to continuously interpolate between the situation where one 

studies a single cascade sample at arbitrary fine scales and where at 
fixed scale, the sample length (number of cascades realizations) be- 
comes infinite. We show that scaling exponents of "mixed" partitions 
functions i.e., the estimator of the cumulant generating function of 



, the cascade generator distribution, depends on some "mixed asymp- 



totic" exponent x respectively above and beyond two critical value 
and pj. We study the convergence properties of partition functions 
in mixed asymtotics regime and establish a central limit theorem. 
These results are shown to remain valid within a general wavelet 
analysis framework. Their interpretation in terms of Besov frontier 
^ ■ are discussed. Moreover, within the mixed asymptotic framework, we 

' establish a "box-counting" multifractal formalism that can be seen 

, as a rigorous formulation of Mandelbrot's negative dimension theory. 

Numerical illustrations of our purpose on specific examples are also 
provided. 

in 
o . 

OO . 1. Introduction. Multifractal processes have been used successfully in 

many applications which involve series with invariance scaling properties. 
Well known examples are fully developed turbulence where such processes 
^ . are used to model the velocity or the dissipation energy fields [6] or finance, 

^ I where they have been shown to reproduce very accurately the major "styl- 

ized facts" of return time-series [U \5\ [21] . Since pioneering works of Mandel- 
brot \14:\ [T5] , Kahane and Peyriere [llj , a lot of mathematical studies have 
been devoted to multiplicative cascades, denoted in sequel as Al-cascades 
(see e.g. refs [HJ [131 [181 [22]). One of the central issues of these studies was to 
understand how the partition function scaling exponents (hereafter denoted 
as To{q)), are related, on one hand, to the cumulant generating function of 
cascade weight distribution and, on the other hand, to the regularity prop- 
erties of cascade samples. Actually, the goal of the multifractal formalism 
is to directly relate the function ro(g) to the so-called singularity spectrum, 
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i.e., the Hausdorff dimension of the set of ah the points corresponding to 
given Holder exponent. Let us mention that recently continuous versions 
of multiplicative cascades have been introduced [U [2] : they share most of 
properties with discrete cascades but do not involve any preferential scale 
ratio and remain invariant under time translation. In these constructions, 
the analog of the integral scale T, i.e., the coarsest scale where the cascade 
iteration begins, is a correlation time. 

In all the above cited references, the main results concern one single cas- 
cade over one integral scale T in the limit of arbitrary small sampling scale. 
However, in many applications (e.g., the above turbulence experiments) 
there is no reason a priori that the length of the experimental series cor- 
responds to one (or few) integral scale(s). From a general point of view, as 
long as modeling a discrete (time or space) series with a cascade process is 
concerned, three scales are involved : (i) the resolution scale / which corre- 
sponds to the sampling period of the series, (ii) the integral (or correlation) 
scale T and (iii) the size L of the whole series. Using these notations, the 
total number of samples of the series is 



Therefore, when modeling a discrete series with a multifractal process, vari- 
ous types of asymptotics for N +00 can be defined. The "high resolution 
asymptotics" considered in the literature, corresponds to Z ^ whereas 
L is fixed. On the other side, one could also consider the "infinite historic 
asymptotics" that corresponds to L — >■ -|-oo whereas I is fixed. If we define 
Nt to be the number of integral scales involved in the series 



Thus, the high resolution asymptotics corresponds to Nt fixed and Ni — > 
-|-oo whereas the infinite historic asymptotics corresponds to Ni fixed and 
Nt +00. But in many applications, it is clear that since the relative 
values of Nt and A'^; can be arbitrary, it is not obvious that one of the two 
mentionned asymptotics can account suitably for situation. This leads us 
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and Ni the number of samples per integral scale 





then we have 



N = NtNi. 
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to consider an asymptotics according to which Nt and Ni go to infinity 
(and therefore goes to infinity) and at the same time preserve their rela- 
tive "velocities", i.e., the ratio of their logarithm. Some of us, have already 
suggested the following "mixed asymptotics" [121 IHl ISO] : 

Nt = N^, 

where x ^ I^"*" is a fixed number that quantifies the relative velocities of Nt 
and Ni. Thus, 

• X = corresponds to the high resolution asymptotics, 

• X ^ +00 corresponds to the infinite historic asymptotics, 

and all other values are truly "mixed" asymptotics. Successful applications 
of the mixed asymptotics have already been performed [19\ [20] . In this paper 
we revisit the standard problems of (i) the estimation of cascade generator 
cumulant generating function in the mixed asymptotic framework and of (ii) 
the multifractal formalism or of how to relate this function to a dimension- 
like quantity. 

The paper is organized as follows: in section [2] we recall basic definitions 
and properties of 7V4-cascades. Section [3] contains the main results of this 
paper. If we define a multifractal measure jl as the concatenation of inde- 
pendent 7V4-cascades of length T, with common generator law W, then we 
show in Theorem [21 

, / JV-l \ 

^^^^ log i^N-' J2 ([^^' + ^P- log2 ■■= r{p) + 1 

for p in some range These critical exponents p^, are related to 

the two solutions, h~, of the equation D(h) = — x where 

D{h) = inf {ph - t{p)} 

is the Legendre transform of r. The convergence rate is studied in Section 
13. 5[ Let us stress that the range of validity on p of this convergence is wider 
in the mixed asymptotic framework (x > 0) than in the high resolution 
asymptotic (x = 0). As a consequence we can relate D{h) to a "box-counting 
dimension" (sometimes referred to as a box dimension [10] or a coarse- 
grain spectrum [23]), and derive, as stated in Theorem [H a "box-counting 
multifractal formalism" for jl 

-^#{k e {0, . . . , iV} I A {[kl, (k + 1)/]) G [l^-^, l>^+^]} ^ /-^W 
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in the range of [h^,h~]. Since for % > 0, D(h) can take negative values in 
previous equation, this can be seen as a rigorous formulation of Mandelbrot's 
negative dimension theory [16^ [TT] . In section [U we extend previous results 
to partition functions relying on some arbitrary wavelet decomposition of 
the process. In Section [5] we give an interpretation of the results connected 
with the Besov frontier associated with our multifractal measure. Finally, 
in section [6] we discuss some specific examples where the law of the cascade 
generator is respectively log-normal, log-Poisson and log-Gamma. For illus- 
tration purpose, we also report, in each case, estimations performed from 
numerical simulations. Auxiliary useful Lemmas are moved to Appendices. 

2. A^-cascades : Definitions and properties. 

2.1. Definition of the Ai-cascades. Let us first introduce some notations. 
Given a j-uplet r = (ri, . . . ,rj), for all strictly positive integer i < j, we 
note r\i the restriction of the j-uplet to its first i components, i.e., 

r\i= {ri,...,ri), Vi € {1, . . . , j}. 

By convention, if j = 0, we consider that r = and in the sequel, we 
denote by rr' the j + j'-uplet obtained by concatenation of r € {0, 1}-' and 
r' € {0, ly . Moreover, we note 

_ ^ |2^' Eii ri2-\ if r / 
^ [0 ifr = ■ 

Let fix T € M"*"* and /c G N. We define Ij^k as the interval 

(2) Ij^k = [k2-^T,{k + l)2-^T]. 

Thus, for any j G N*, the interval [0, T] can be decomposed as 2^ dyadic 
intervals : 

[o,r]= U Ij^r. 
re{o,ip 

Let us now build the so called 7V4-cascade measures introduced by Mandel- 
brot in 1974 [15]. Let {I^r}rG{o,i}J, jen* ^ of i-i-d random variables of 
mean E [Wr] = 1. Given j G N*, we define the random measure fxj on [0, T] 
such that, for all r G {0, 1}-', the Radon-Nikodym derivative with respect to 
the Lebesgue measure is constant on Ij^r with: 

(3) ^ = f[Wr\^, on4^,forrG{0,iy. 

i=l 
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As it is well known [TT], the measures /Xj have a non-trivial limit measure 
/ioo, when j goes to oo, as soon as E[T4^1og2VF] < 1. Moreover, the total 
mass 

i 

//oo([o,T]) = hm r2-^^ J2 UWri,, 

rG{0,lp «=1 

verifies E [//qo ([Oj^])] = T. Let us remark that if r € {0, 1}-' then by con- 
struction we have: 




where is a 7V4-cascade measure on [0, T] based on the random variables 
Wrr' for r' G Uj>i{0, !}•' . This equality is usually referred to as "Mandelbrot 
star equation". 

In the sequel we need the following set of assumptions: 



(5) E[Wlog2W]<l, F{W = 1) <1, 

(6) ¥{W > 0) = 1, E [WP] < 00 for aU p£R. 

Let t{p) be the smooth and concave function defined on M by 

(7) T{p)=p-log^E[WP]-l. 



Let us notice that log2 E [W^] is nothing but the cumulant generating func- 
tion (log-Laplace transform) of the logarithm of cascade generator distribu- 
tion. It is shown in [11] that for p > 1, the condition r(p) > implies the 
finiteness of E [fi^o ([0,^])^]. By Theorem 4 in [18], the conditions ^ imply 
the existence of finite negative moments E [fiQ^ ([0, T])^], for all p < 0. 

2.2. Multifractal properties of A4- cascades. A M-cascade is a multifrac- 
tal measure and the study of its multifractal properties reduces to the study 
of the partition function 

(8) s,{j,p) = J2 ^oo(/i,fcr. 

k=0 

Basically, one can show \18\ [22] that, for fixed p, this partition function 
behaves, when j goes to 00, as a power law function of the scale | = T2~^ . 
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More precisely, let us introduce the two following critical exponents: 

Pq = inf{p > 1 I pt'{p) — t{p) < 0} € (1, oo] 
Pq = sup{p < I pt'{p) — t{p) < 0} G [— oo, 0). 

If (resp. Pq ) is finite we set = t'{pq) (resp. = t'(pq )). 

Theorem 1. Scaling of the partition function 12^ 

Let p G M, the power law scaling exponent of Sfj_{j,p) is given by 

(9) lim f^, >tq{p), 

j^ca —J as. 

where To(p) is defined by 

( t{p), Vp G (Po ,Po ) 

(10) roip) = l h+p, \fp>p^_ 

[ K ^ Po 

The proof can be found in [22]. The convergence in probability of ([9]) was 
obtained in the earlier work [18]. This theorem basically states that S^{j,p) 
behaves like 

5^(j,p)=.2-^--"(p). 

Let us note that the partition function ^ can be rewritten in the following 
way 

(11) Sf,{j,p)= ^^oo{I,,rY 

rG{0,lp 

and using one gets 

(12) 5^(i,p) = 2-^> ri^.l^(M£^([0,r]))', 

re{0,l}J «=i 

where the {/Z^^([0, r])}^g |o,i}j are i.i.d. random variables with the same law 
as /ioo([0,T]). Thus, a simple computation shows that 

E [5^(j,p)] = 2-^>2-'E [W^y E [/t^oo([0, T]Y] = 2-^<p^¥. [/.oc([0, r])^ . 

One sees that the last theorem states that, in the case p G [P(^,Pq'], S^{j,p) 
scales as its mean value. 

On the other hand, the fact that for p ^ [Pq ,Po^] the partition function 
scales as given in (jlOp instead of scaling as its mean value, is referred to 
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as the 'linearization effect'. A possible explanation of this effect is that for 
p larger than the critical exponents Pq (resp. smaller than Pq), the sum 
involved in the partition function ([TT]) is dominated by its supremum (resp. 
infimum) term. Thus one should not expect a law of large number to hold 
for the behavior of this sum. Another possible interpretation of this theorem 
in the case p > Pq is given in Section [5l 

3. Mixed asymptotics for Al-cascades. 

3.1. Mixed asymptotics : definitions and notations. A convenient way to 
construct a multifractal measure on M_|_, with an integral scale equal to T, is 
to patch independent realizations of 7W-cascades measures. More precisely, 
consider {/i^^jmGN a sequence of i.i.d A^-cascades on [0,T] as defined in 
Section [2.11 and define the stochastic measure on [0,oo) by: 

(13) /i ([ti, t2]) = E /"^'^ - *2 - mT]) , for all < ti < ta- 

m=0 

This model is entirely defined as soon as both T and the law of W are fixed. 
The discretized time model for the N samples of the series is {/u[A;/, (k + 

l)^]}o<A:<Af-l- 

3.2. Scaling properties. In this section, we study the partition function 
for the measure /x as defined in Eq. (jl3p in the mixed asymptotic limit. T is 
fixed, we choose the sampling step 

/ = T2-^, Ni = 2^ 

and the number of integral scales is related with the sampling step as 

Nt = [Nl^\ ~ 2^^, 

with X > fixed. According to ([T]), one gets for the total number of data: 

N = Nt2^ ~ 2J'(^+x). 

The mixed asymptotics corresponds to the limit j +oo. The partition 
function of /i can be written as (recall Q): 

(14) Sf,{j,p) = Y.~^ihkY 

k=0 

(15) = E s!r\j,p), 

m=0 
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where S^\j,p) is the partition function of i.e., 

(16) s^;^Hj,p) = ''£f,tHkkr- 

k=0 

Let us state the results of this section. We introduce the two critical expo- 
nents in the mixed asymptotic framework: 

(17) p+ = inf{p > 1 I pt'{p) - t{p) < -x} S (1, oo] 

(18) p~ = sup{p < I pt'{p) - t{p) < -x} G [-00, 0), 

and we set when these critical exponents are finite /ij = t'(pJ), h~ = t'{p~). 

Theorem 2. Scaling of the partition function in a mixed asymp- 
totics 

Let p G M and fi be the random measure defined by (I13p where the law of 
W satisfies ([SI)-®. We assume that, either p'^ = oo, or p+ < oo with 
r(pj) > 0. Then, the power law scaling of Sfi{j,p) is given by 

(19) limi^iiMlZ)_,^(^), 
where T^{p) is defined by 

{r{p)-X, ^P^{Px^Px) 
h+p, '^p>p^ 
h~p, Vp < p~ 

Remark 1. If p'^ = oo then simple considerations on the concave func- 
tion T shows that t[p) > for all p > 1, and hence the cascade measure has 
finite moments of any positive orders. Otherwise the assumption t(pJ) > 
is stated in Theorem\^to insure E [^oo([0, T])*'] < oo for p € [0,pj). Such 
assumption was not needed in TheoremUl since on can check that necessarily 
r(p+) > 0. 

Remark 2. Let us stress that the behavior of the partition function is 
largely affected by the choice of a mixed asymptotic: the 'linearization effect' 
now occurs for p in the set {—oo,p~) U (pJ,oo), which is smaller when x 
increases. 
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Theorem 3. Scaling of the supremum and the infimum of the 
mass in a mixed asymptotics 

Assume ([5]) -([6]). Then, if < oo, one has, 

(20) lim — — — ^] -^'^-^ _ /j+ almost surely, 
and if > —oo, one has, 

(21) lim — — ^'^[0'^ ^1 ^( J'^-^ _ ^ almost surely. 

j^+oo — j ^ 

The theorem [3] shows that when the 'linearization effect' occurs, the scal- 
ing of the partition function (jl4p is governed by its supremum and infimum 
terms for respectively large positive and negative p values. 

These theorems will be proved in three parts. In Section 13.3.21 we will 
prove Eq. (fT9]) of Theorem [2] only for p € {p~,p^). In Section [3.3.31 we will 
prove the case p ^ (p~,p^) and Theorem [3] is shown in Section [3.3.41 to be 
a simple corrolary of this last case. 

3.3. Proof of Theorem\^ First we need an auxiliary result which is help- 
ful in the sequel. 



3.3.1. Limit theorem for a rescaled cascade. For each m we denote as 
(VKj'"^)^gUj{o,i}J the set of i.i.d. random variables used for the construction 

of the measure Moreover we assume that for each m > 0, j > 0, 

r € {0, 1}-^ we are given a random variable Z^"^'^\ measurable with respect 
to the sigma-field a (wj:!^^ \ r' G Uj{0, ly^. We make the assumption that 

the law of Z^"^'^^ does not depend on (m, r), and denote by Z a variable 
with this law. 

Let us consider the quantities, for p G M: 

(22) Mf'\p) = 2~^P J2 f[{w^\^Y ^^"^^'^^ 

rG{0,lp '=1 

and 

Nt-I 

(23) AfM= E M^rHp)- 

m=0 
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Proposition 1. Assume that for some e > 0, E \Z 
—Pt'{p) +t{p) < X, then: 

2i(^(p)-x)_Yj(p) E [Z] , almost surely. 

Proof. From (l22D-([23]) and the definition ([7]) we get, 
E [Mj{p)] = Nt2^2-^pK [WPy E [Z] 

-j^oo m-^^^p^E [Z] . 

Hence the proposition will be proved if we show: 
(24) 2^(r{p)-x) (_^. (p) _ E [_^^. (p)] ) 



< oo anc 



0, almost surely. 



For an arbitrary small e > 0, we study the L^^'^(F) norm of the difference. 
Set, 



(25) 



\Mj{p)-E[M,{p)] 



L]^' < C2^^E 



Applying successively lemmas [T] and [2] of Appendix [Aj we get: 

\Mf\pf' 
j 

< C'2~-''[(^+')^(P)~^] 2-*^^(p(i+^))2^(^+^)^(P^ 

A:=0 

We deduce that 2J(^(p)-^)(^+^)Ly"' is bounded by the quantity: 

j 

fc=o 

Clearly, as soon as 2-x^2-^(p(^+^))2(i+')^(p) < 1, this quantity is, in turn, 
bounded by C2~^^^ for some e' > 0. Taking the log, a sufficient condition is 

r{p){l + e)-T{p{l + e)) ^ 

< X 

e 

which is implied for e small enough by 

-pr'ip) +t{p) < X- 

Thus we have shown that 2^^'^^^^"^^^^'^''^ L]^"" is asymptotically smaller than 
2~^^ with some e' > 0. Using the Bienayme-Chebyshev inequality leads to 



!){2i(-(p)-x)|AAj.(p) _ E [Mj{p)] \>v}< 



for any 77 > 0. A simple use of the Borel Cantelli lemma shows 



< 



□ 



imsart-aap ver. 2007/12/10 file: paper_ver5-3_IMS.tex date: May 2, 2008 



11 



3.3.2. Proof of Theorem IE for p G {p^,Px)- From (fT5]) - (fT6]) and the 
representation (jl2p for the partition function of a single cascade, we see 
that that Sfi{j,p) exactly has the same structure as the quantity Mj{p) of 

section [3.3.11 where Z^"^'"^^ = 7I^'^^([0, T])^ are random variables distributed 
as Z = /ioo([0,T])P. 

By definition (recall (fT71) - (fT8]) ). the condition —pr'ij)) + r(p) < x holds 
for any p S {p^,P^), and by Remark[Tl E [iZl^+^j < oo for e small enough. 

Thus, an application of Proposition [1] yields the almost sure convergence: 

(26) 2^--x(p)5-(j-p) = 2^-«rt-x)5-(j-,p) i^E[/ioo([0,r]r]. 
This proves the theorem [2] for the case p € {p^,p^)- 

3.3.3. Proof of Theorem for p ^ {p~,p^). The following proof is an 
adaptation of the corresponding proof in [23]. We need the following nota- 
tions: 

Sf*{j)= sup fioo{[k2-^T,{k + l)2-^T]), 

k£[0,N-l] 

m..,ip) = limsup ]2MlMA, ^^„,(,) = n^inf ]2^lMA, 

j^oo -3 -J 

log S-(j)* log S-(j)* 
m* = hm sup — - — , m*„ f = hm inf — - — . 

j^oo -3 i^°o -J 

In Section 13.3.21 we proved that for all p S {p^,P^) the following holds 
almost surely: 

nisupip) = niinfip) = Tx(p)- 

We may assume that on a event of probability one, this equality holds for 
all p in a countable and dense subset of (p~,p+). 
From the sub-additivity of x i-^ x'', 

vp e]o, i[, vp G M s^ipjy < Sf^ippj), 

and thus 

, s . ruinfipp) 

rriinfiP) > • 

P 

But we have seen that minf{pp) = T^{pp), for a dense subset of pp G {p~,p^). 
Assume now for simplicity that p > p'^ and let p {p^/p)^ S^t 

(27) Vp>p,^ !!!hizM>!xM)^iM^ = /,j, 

P Px Px 
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where we have used (1171) 

On the other hand, let p > 0, g € [0,pj), and (f G [0, g), we have 

Sf,{j,q) = J2 f'oci[k2-^T,{k + l)2-^T]y 

k=0 



q-q 



Thus 



then 



rusupiq) >{q- q') ^"''^^^^ + rriinfiq'), 



msupjp) ^ rrisupiq) " ruinfiq') ^ Txil) ~ '^x(gO 
p ~ q — q' q — q' 



Taking the hmit q' ^ q 

P ge[o,p+) 
Merging this last relation with (j27p leads to 

(28) Vp>p+, ^+ < < < , 

which proves Theorem [2] for p € [pj, +oo[. The proof for p < is similar. 

3.3.4. Proof of Theorem\^ The following proof is an adaptation of the 
corresponding proof in [23]. We have for p > 0, 

Sp*{jY < S^{j,p) < NSp * Ur = [2^^\2^Sf,*ijr, 

thus 

P^T'inf,sup > mnf,sup{p) > 1 + X+Pm*nf,supi 

which means that 

^ ^* ^ nT'inf,sup ip) 

— ^inf,sup — ! 
p p ^ p 

and taking the limit p — > +00 and using (j28|) proves that 

* * J + 

''^sup ~ ''^inf ~ '^x 1 

which proves ([20]) . The proof of ([2T]) is obtained analogously by considering 
p < 0. 
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3.4. Multifractal formalism and "negative dimensions". Let be the 
Legendre transform of r(p) : 

D(h) = miniph — T(p)), 
p 

The multifractal formalism [7] gives an interesting interpretation of D{h), 
as soon as D{h) > 0, in terms of dimension of set of points with the same 
regularity. For A^-cascades, this formahsm holds [H], i.e., D(h) corresponds 
to the Hausdorff dimension of the points t € [0, T] around which fi^o scales 
with the exponent h : 

(29) D{h) = d^mn U hmsup ^^^2 l^oo{[t - e t + e]) ^ ^ 

I e^O log2(e) J 

The r.h.s. of (j29p is usually referred to as the singularity spectrum and there- 
fore the multifractal formalism simply states that D[h) can be identified with 
the singularity spectrum of the cascade. 

In a mixed asymptotic framework, our next result shows that some kind 
of multifractal formalism still holds for D{h) < in the sense that D{h) 
governs the behavior of the population histogram per sample of measure 
values at scale as estimated over 2^^^ cascade samples. In other words, 
D[h) coincides with a box-counting dimension (sometimes referred to as a 
box dimension [10] or a coarse-grain spectrum [M])- Hence the Legendre 
transform of t{p) can be interpreted as a "population" dimension even for 
singularity values above and below and . Since for these values, one has 
D[h) < they have been called "negative dimensions" by Mandelbrot [16] . 
This simply means that they cannot be observed on a single cascade sample 
but one needs at least 2^^ realizations to observe them with a "cardinality" 
like 2^^^~^^^'^^\ In that respect, they have also been referred to as "latent" 
singularities [T7]. 

Theorem 4. Assume p'^ < oo, p~ > —oo and t{p^) > 0. Let h G 
{h+,h~), then: 

(30) limlim-log# {yt G {0, . . . ,iV - 1} ! 2--'('^+^) < fl{Ij^k) < 2'^^^-^^} 

= X + D{h), 

(31) limIE^-log# {A; G {0, . . . , iV - 1} [ 2~J('^+") < fi{Ij^k) < 2~^'^''~'^} 

= X + D{h). 
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Proof. 1^* step: We focus on the cases that yield to negative dimen- 
sions, i.e. h € {h^,h^) U (/iq For simphcity assume h £ = 
{t'{p^),t'{pq)). We can write h = t'{p) for some p G {Pq ^P'^) and if we 
define 

= -p-r'ip) G (0,x) 

we easily get that h = t'(pJ) and D[h) = —x'- Thus the theorem amounts 
to assess the magnitude of 

# {a; € {0, . . . ,iV - 1} I 2-^'("'(Px')+^) < ^(/^.,) < 2-^'(^'(^x')-^)| 
as ~ 2^'-^-^'\ 

First we derive a lower bound for this cardinality. The idea is to split the 
data into blocks of size 2-' 2^^' and rely on the behavior of the supremum 
of fl{Ij^k) under mixed asymptotic with index x' ■ More precisely let N' = 



2^ 



and define the blocks 
B„ = {aN\ . . . , (a + l)N' - 1}, for a = 0, . . . , M - 1 := [N/N'\ - 1. 
Fix a € {0, . . . ,M — 1}, then for any pi < P2 < pt' we have 

(32) sup Khkr-'' > g^^^" ^l^^-'^i"' := 2^(v(^^0-v(p^))q('^)(pi,p2), 



where 



Clearly the law of Q^j^\pi,P2) does not depend on a, andgf(pi ,P2) is the 
ratio of two rescaled partition functions in mixed asymptotic with index x' ■ 
Hence by ([26|) . Qf'\pi,P2) converges almost surely to the non zero constant 
E [fioo{[0, T]y^] /E [fioc{[0, T])Pi]. We deduce that for all fixed a < M-1, and 
P1tP2 £ (0,pj/), the sequence (l/<3j"^(pi,P2))j>o is bounded in probability. 
Write now, by ([32]) . and r^'ip) = t{p) — x' for p < p^,, 

sup Kk,)2^^''^''P^'^ > Qf\p^,P2)2-^^''^-^'^<'^-'^ 

fcGB„ 

and choose pi, p2 fixed but close enough to p^/. This yields, for all j > 0: 
sup f,{I,,,)2'^^''-^>'^ > Qf{p,^p2Wl\ 

fcGB„ 
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Then using that 1 / Q^"^ (pi , P2) is bounded in probabihty, we get P(sup;jgiB^ fj'{Ij,k) > 

> F{l/Qf\pi,p2) < 2-''^/2) > 1/2 for j large enough. 
Remark now that the cardinality of the set {A; € {0, . . . , A^ — 1} | jl{Ij^k) ^ 
2 ^ ^ x" '} is immediately lower bounded by the sum 

M-l 



El 



a=0 S supfegB u(/j^fc)>2 X 



of i.i.d. Bernoulh variables with parameter greater than 1/2. Then it is easily 
deduced, using the Borel Cantelli lemma and M ~j^oo 

2(x~x')j that with 

probability one: 

(33) lim2~(>^"'^')-'# jfc G {0, . . . ,iV - 1} I fi{Ij,k) > 2"^'^"'^^x')+") | > 1/4. 

We now focus on upper bounds for the cardinality of the set 
{ke{0,...,N-l}\ > 2-^'("'(^x')+^)} 



where tj is some real number in a neighborhood of zero. It is simply derived 
from the connection with the partition function that for any p > this 

cardinahty is lower than »S^(j,p)2"'^^ Applying this with p = p^, 

and since, by ([IT]), t^{p^,) = - X = pr'{p'^>) + x' - X we get the 

following upper bound, 

(34) #{ke{0,...,N-l}\ Ji{I,^k) > 2-^("'(<'^+'')} 

< 5^ (j, p+ ) 2^"^ K' ) 2^^^-^'^<' ) . 

By p^, < p+, the convergence result (j26]) with p = p^, apphes and we deduce 
for r] = —€ < that: 

(35) Iim2-(>^-^')-'#{A; G {0, . . . , iV - 1} | fl{Ij^k) > 2~-''^"'^^x'^"'^} = 0. 

Then ([30|) is a consequence of (p3|) and ([35]) . 

Finally, the upper bound ([^T]) is directly obtained by applying ([M]) with 
= e. 

2"*^ step: We now deal with the more classical case h G [/iq',/iq]. It is 
known, from the multifractal formalism for a single M-cascade on [0, T] (see 
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[3l[24]), that with probability 1: 

(36) lim lim-log#{A: g {0, . . . , 2^' - 1} | G [2-i(''+-), 2-^(''-^)]} 

= D{h). 

For m G {0, . . . , A^t — 1}, and e > 0, r/ > 0, denote by j4^"^^(7/, e) the event: 
# {fc G {0, . . . , 2^- - 1} I Klj,2.^+k) e [2-^-(^+-), 2-^-(^--)]} > 2^-(^W-''). 



Using the independence of the A4-cascades {fj,^^^)m., these events are inde^ 



pendent and P(y4[™^ (77, e)) does not depend on m. Moreover, by ([36]) . for 



any r/ > 0, there exists e > such that F{A^j^\ri,£)) - — ^ 1. We easily 
deduce that for any 77 > and e small enough: 

lim-— If^C"")^ M > 1/2, almost surely. 

J ^ m=0 

Since the cardinality of {/c G {0, . . . , - 1} | /i(/,-,fc) G [2-J('^+^), 2-J('^-^)]} 
is greater than 2-'(^('*)~'') X]m=o^ ^{^''"'(r; e)} deduce that the left hand 
side of (jSOp is greater than D{h) +x~ "Hi for any r/ > 0. To end the proof, it 
suffices to show that the left hand side of (|3ip is lower than D{h) + x- This 
is easily done, as in the end of the first step, by relying on the asymptotic 
behavior of the partition function. □ 



3.5. Central Limit Theorems. In this section, we briefly study the rate 
of the convergence of Sfi{j,p) as j — > 00 of (fT9]) in Theorem [2j Using the 
same notations as in the proof of Theorem [2l we write: 

(37) S^{j,p) = [2^^\ 2~^^(P^K[f,^{[0,T]r]+A^+Bj 

where 

Nt-I ( j \ 

^i=E2"^'' E n^.^^'^n (^"'^''H[0'^])'-iE[^oo([o,r])^]) 

m=0 re{0,l}J \«=1 / 



and 

^ fo~jp \ " TT w 

^1 



^j=E(2"^'" E flw^/?^-2-^'^(p))E[^oo([o,r])^ 



m=0 rg{0,l}J «=1 
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Proposition 2. Assume ([S])-® and that, either p'^ = 00, or p+ < 00 
with t{p+) > 0. If < p < _p+/2 then, 

2iW2p)-x)/2^^. _^(0, Var(^oo([0,T])P)). 

Proof. Consistently with the notations of Section [3.3.H we define, for 
every r G {0, Ip and m = 0, . . . , A't^ — 1, the random variables Z^*"'^-* = 
/i^''^([0, r])P-E [/ioo([0, T])P] and denote by Z = ^ioo([0, r])P-E [;Uoo([0, T])p] 
their common law. Furthermore, we will need the quantity 

(38) yim,rM = 2^(-(2rt-x)/22-.P ^^Mf^ ^(-,0^ 

and the following family of fj-fields: for j > 

:= |r| < i, m = 0, . . . , A^T - l) 

and for every A; = 0, . . . , n(j) = 2^[Nt — 1) 

^fc^^. = V (j(Z(™'^), f + 2Jm < A:). 

For fixed j, we have a one-to-one correspondence between {m,r) and A; = 
f-f 2-^m, so abusing notation slightly, we write rj^jip) instead of r]m,r,jip) in 
psp when no confusion is possible. With these notations, 

k=0 

where rj^jip) is ,,-measurable and 

^ [Vk,j{p) I ^k-i,j] = 0, V/c = 0, . . . , n(j). 

Thus, we are dealing with a triangular array of martingale increments. Let 
us consider the sum of the conditional variances: 

(39) V, = Y,E[r^k,j(j^f\:Fk-i,, . 

fc=0 

We have 

Nt-1 j 

V, = Var(Z)2^(-(2p)-x)2-2.P J] ^ J] (^r^O ' ' 

m=0 rG{0,l}J «=1 
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thus by application of Proposition [T] (with the choice of Z^"^'^^ equal to 1) 
we get, 

Var(Z). 

Hence the proposition will be proved, if we can show that the triangular 
array satisfies a Lindeberg condition: for some e > 0, 



J-+00 



A:=0 

But, we have 



^ 0. 



re{0,l}J «=i 



vr = E 



\w 



2+e 



and by application of the Proposition [H the order of magnitude of vj^^ is 
2j((r(2p)-x)(i+^/2)-(r((2+p)e)-x)) ^ ^hus, it Can be seeu that v}'^ converges to 
zero, for e small enough, by the condition 2pT'{2p) — T{2p) > — x- This ends 
the proof of the proposition. □ 

Proposition 3. Assume ([S])-® and that, either p^j^ = oo, orp+ < cx) 
with t(p+) > 0. Then: 

1. IfT{2p) - 2t{p) > we have, 

(40) 2^{r{p)~xl2);Q. ^ _^(o, C(p)) , 

where c{p) > depends on the law of W and p. 

2. IfT{2p) - 2t{p) = 0, we have Var{Bj) = 0(j2-j(^(2p)-x))_ 

3. IfT{2p) - 2t{p) < 0, we have Var{B-) = 0(2-^(^(2p)-x)), 

Proof. Denote z^j"*^ the measures defined at the step j of the construc- 
tion of the M-cascade on [0, 1] based on W^^/E \WP]: 



'f)([0,l])=2-^- n(VF(;f''^fE[VFP]-^ forme{0,...,iVT-l}. 

re{0,lp i=l 



With this notation we have, 

Nt~1 

(41) i?,- = E[^oo([o,r])^]2-^'-(^) (^r^([o,i])-i) 

m=0 
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and using Kahane and Peyriere results [TT], we know that for each m the 
sequence {i'j"^\[0, is bounded in L'^ as soon as T{pq) — qr^p) > 0. 
We first focus on the case t{2p)—2t{p) > 0. Hence the sequence (^'j"*^([0, 

is bounded in L^'^'^-norm for some e > 0. Using that (z/j™^([0, 1]) — l)m is a 
centered i.i.d. sequence and classical considerations for triangular array of 
martingale increments, one can show that a central limit theorem holds: 

N~'/' Y: (-r^([0,T])-l)^Ar(0,Var(.W([o,l]))), 

m=0 

where u^\[0,l]) = liuij^oo From this and (jH]) we deduce (00]) 

with c(p) =E[^ioo([0,r])P]2Var(z.S'([0,l])). 

In the cases r(2p) — 2r(p) < 0, by (0T]) again we have 

Var(i?,-) = 2-^'2-(p)E[^oo([0,r]f]2Var(z.f ([0,1])). 

Now Var(z^j°^([0, 1])) = E k-°^([0, l])^] - 1 is unbounded as j oo, but a 
careful look at the computations in Lemma [2] with e = 1 yields to 



E 



z.f([0,l])2 ~,^ooE2"'2'(2^^^^-^ 



1=0 



We deduce that Var(Sj) = 0(2-J(2^(p)-^) 2-'2'(2^(p)-^(2p))-)_ r^^^^^ ^j^g 
theorem follows in the cases r(2p) — 2r(p) = and T{2p) — 2r(p) < 0. □ 

Remark 3. By ([37|) t/ie difference between 2^'^^p^~^'^^ S^{j,p) and its 
limit is decomposed into two dissimilar error terms: particularly the fact 
that the contribution of Bj converges to zero is due to the observation of a 
large number of integral scales, whereas the contribution of Aj vanishes as 
the sampling step tends to zero. 

In the case r(2p) — 2r(p) > 0, the contribution of Bj strictly dominates 

and 2-(^(P)->^)5^(j,p) -E[/Xoo([0,T])P] is of magnitude 2"^'^ ~ iV^^^^^ 

If t{2p) — 2t{p) < 0, the magnitude of Aj and Bj are the same and 
2-('^(p)~x)5-(j^p) — E [/Xoo([0, T])P] is asymptotically bounded by terms of 
magnitude 2J'/2(-x+2T(2p)-r(2p)) ^ This rate of convergence is slower than Nj,^^'^ . 

4. Extension to wavelet based partition functions. The behavior 
of the partition function provides an evaluation for the regularity of the 
sample path of the process t /i([0, t]). However, it is more natural to 
assess this regularity via the behavior of wavelet coefficients. 
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4.1. Notations. In this section we assume, for notational convenience, 
that T = 1. Consider now g a "generaUzed box" function. It is a real valued 
function that satisfies the following assumptions 

(HI) g has compact support included in [0, 2'^], for some J > 0. 

(H2) g is piecewise continuous. 

(H3) g is at least non zero on an interval. 

Following the common wavelet notation, we define 

gj^kit) = g{2H - k). 

The support of gj^kit) is 

(42) Supp gj^k = i^-'k, 2~^k + 2^"^]. 

In the sequel, if /i is a random measure, for any Borel function / we will use 
the notation 

(/.,/)= / f{t)d^l{t). 



4.2. The generalized partition function : scaling properties. We define 
the generalized partition function of an M-cascade /Xqo on [0, 1] at scale 2~^ 
as 

21-2^-1 

(43) S,,g{j,p)= J2 K/^oo,<7,,fc)r. 

fc=o 

Remark that for simplicity we removed a finite number of border terms, and 
that, in the case g{t) is the "box" function g{t) = l[o,i](i) we recover the 
partition function of Section 12.11 

Let us study the scaling of E [Sfj,^g{j,p)]. 

Proposition 4. Assume ([S])-®. Then, we have Ki2-^'^'^p'^ < E [5^,g(j,p)] < 
K22~^'^^P^ for Ki, K2, two positive constants depending on p, W and g. 

Proof. Since \g{t)\ is clearly a bounded function, we have 

E[|(^oc,5i,fc)n < CE [^U[2-'k,2-^k + 2J-^]r 

where C is a constant. We write /Uoo([2~-'/c, 2~-^/c+2"^~-^]) = ^21=0^ k''oo{Ij,k+i) 
and deduce 



(44) E[|(^oo,ffj,fc)n <CE |/ioo[0,2 



-J IIP 
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where K only depends on g and the law of W . By (|43l) we get the upper 
bound for E [S^^g{i,p)]. 

For the lower bound, let us write that S^^g{j,p) is greater than 

X! i(/^oo,5j,2-^fc')r- 

fc'=0 

But gj 2Jk' is supported on [A;'2'^~-', + 1)2"^~-'], thus applying Lemma[3]in 
the Appendix [B] with a = j — J, we deduce: 

where, in law, Z is equal to {fJ-oo, gj,o)- Thus E [|(/^ooi <?j,fc)|^] is greater than 
2p(J-^)e [WPy-^ E [| (/.oo, gj,o) r] = K2-^-(-(^')+^)E [| (/ioo, 5J,o) T] • 

Applying Lemma[4]with / = gj^ shows that E [\{Hoo, gj,o)\^] is some positive 
constant. Then the lower bound for E [S^^g{j,p)] easily follows. □ 

4.3. The partition function in the mixed asymptotic framework. Follow- 
ing (lisp , we define the partition function in the mixed asymptotic framework 
as 

Nt~1 
m=0 

where sj^\j,p) is the partition function of i.e.. 



2J_2'^-l 

E 

fc=0 



We have the following result. 



Theorem 5. Scaling of the generalized partition function in a 
mixed asymptotic 

Let p > 0, then under the same assumptions as Theorem the power law 
scaling of Sjx,g{j,p) is given by 

1-^ log2WM^ 
j— >oo —j a.s. 
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Proof. Using Proposition|4]we have, hm^^^oo log2 E [Sji^g{j,p)] = t^{p), 
and we just need to prove that almost surely, 

5^,3(j,p)-E[5^,,(j,p)]=o(2-^--x(p)). 

Using lemma [5] and El of appendix [HI this is done in the exact same way as 
the proof of ()24p in Proposition [TJ □ 

5. Link with Besov spaces. Following [9], one may define for a mea- 
sure /i on [0,T], the boundary of its Besov domain as the function : 
(0, oo) — > R U {oo} given by 



2^-1 



i/p 



s^{l/p) = sup ^ (7 G M I sup2^'^ J2 l^(^i,fc)r < ^ 



j>0 



k=0 



The following proposition can be shown (see [9]). 

Proposition 5. The function is an increasing, concave function, 
with a derivative bounded by 1. 

Let us stress that the condition s'^{l/p) < 1 is a simple consequence of 
the Sobolev embedding for Besov spaces. The Theorem [1] characterizes the 
Besov domain for a M-cascade on [0, T]: 



Vp>0, s^^{l/p) 



t{p)+i if i > J_ 

P P po 

h+ + l if i < X 

u ' p p — po 



If we denote = "^^^p"*"^ then, it is simply checked that the condition 
1/p > ^/Pq is equivalent to s'{\/p) < 1. Hence Proposition explains why 
for 1/p < 1/pq the boundary of the Besov domain must be linear with a 
slope equal to one. 

In mixed asymptotic the support of the measure grows with j but we can 
still define, using the notations of Section O the index: 

sf{l/p) = sup j G E I sup2J'^ iN^^2-^ J2 Wj,k)\n < oo 

Then, it is simply checked that Theorem [2] implies s^{l/p) = s{l/p) when 
s'{l/p) < 1 + X) and s^{l/p) = /ij + otherwise. This shows how the 
linear part in is shifted to larger values of p under the mixed asymptotic 
framework. 
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500 1000 



Fig 1. Three synthetic samples of M-cascades with T — 2^'' and — 0.2; (a) log-Normal 
sample, (b) log-Poisson sample with S = —0.1 and (c) log-Gamma sample with (3 — 10. In 
fact we used /ij,„ax+5[w, n + l]„=o l el's a proxy of fioc [n,n + 1] with jmax = log2(r) = 13 
(see Eq. ([3]) for the definition of fij). 
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6. Numerical examples and applications. Our goal in this section 
is not to focus on statistical issues and notably on precise estimates of mul- 
tifractal exponents from empirical data. We rather aim at illustrating the 
results of theorem [2] on precise examples, namely random cascades with re- 
spectively log-normal, log-Poisson and log-Gamma statistics. For the sake 
of sirnplicity we will consider exclusively scaling of partition function for 
P > Qj- In order to facilitate the comparison of the three models, will 
represent the so-called intermittency coefficient, i.e.. 



where t{p) is defined in Eq. ([Tj). This value will be fixed for the three con- 
sidered models. Let {Wr} be the cascade random generators as defined in 
Eq. §i) and let = In Wr- 

In the simplest, log-Normal case the {w^jr are normally distributed ran- 
dom variables of variance ln(2). Thanks to the condition E [Wr] = E [e"^-] = 
1, their mean is necessarily — A^ln(2)/2. In that case, the cumulant gener- 
ating function t{p) defined in Eq. ([7]) is simply a parabola: 



In the log-Poisson case, the variables ujr are written as ojr = mQ ln{2)+5nr 
where the Ur are integers distributed according to a Poisson law of mean 
7ln(2). It results that t{p) = p{l — mo) +7(1 — e^^) — 1. If one sets r(l) = 
and t"(0) = — A^, one finally gets the expression of t(j>) of a log-Poisson 
cascade with intermittency coefficient A^: 



In third case the variables cOr are drawn from a Gamma distribution. If x is 
a random variable of pdf /3"''^^^^x'^''^(^)~^e~^^/r(a ln(2)), then one chooses 
Ur = X + iTiQ ln(2) and it is easy to show that t{p) is defined only for p < (3 
and in this case r(p) = p(l — mo) + a(l — p/P)- By fixing r(l) = 1 and 
r"(0) = A^, on obtains: 



Notice that one recovers the log-normal case from both log-Poisson and 
log-Gamma statistics in the limits 6^0 and /3 +oo respectively. 

Numerical methods for estimating r(p) for p < are trickier to handle 
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For the 3 cases, one can explicitely compute all the mixed asymptotic 
exponents as functions of x: In particular the values of read: 



Pt.a = P 



1 I ^ + X l+W(±,-e ^) 

^ A2/52 



where suffixes n,p,g stand for respectively log-normal, log-Poisson and log- 
Gamma cascades and W{±,z) represent the two branches of the Lambert 
W{z) function, solution of W{z)e^^^^ = z that take (respectively positive 
and negative) real values for the considered arguments. For log-Poisson and 
log-Gamma cases, we have also indicated the asymptotic behavior in the 
limits 5 —>■ and P ^ oo. The values can be easily deduced form their 
definition: = t'{p^). 

In Fig. [1] is plotted a sample of each of the three examples of A^-cascades. 
We chose T = 2^^ and = 0.2 for all models while, in the log-Poisson 
case we have set 5 = —0.1 and /? = 10 in the log-Gamma model. In each 
case, an approximation of the A^-cascade sample is generated. We chose to 
generate ^is (as defined by Eq. ([3])) so that the smallest scale involved is 
f'min = 2~^^T = (we have checked that the results reported below do 
not depend on Imin)- ^n approximation of jl is generated by concatenating 
i.i.d. realizations of /^ig. Then, for each model and for each chosen value 
oi X, Txip) (p = . . . 6) was obtained from a least square fit of the curve 
log2 Sji{j,p) versus j over the range j = . . . 6. Let us recall that, for each 
value of j, the mixed asymptotic regime corresponds to sampling p, at scale 
I = 2~^T and over an interval of size L = 2^^T. The exponents reported in 
figs. [2] and [3] represent the mean values of exponents estimated in that way 
using N = 130 experiments. 

The log-Normal mixed asymptotic scaling exponents for x = 0, 0.5, 1 are 
represented in Fig. [2j For illustration purpose we have plotted T-^{p) + X 
as a function of p: one clearly observes that, as the value of x increases, 
the value of p^ below which the function is linear, also increases while the 
value of the slope /i+ decreases. As expected, when x increases T^{p) + X 
matches t{p) over an increasing range oi p values. Notice that the estimated 
exponents are very close the analytical predictions as represented by the 
dashed lines. Error bars on the mean value estimates are simply computed 
from the estimated r.m.s. over the 130 trials and are reported only for the 
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Fig 2. Estimates of the function r-^{p) of log-normal cascades with = 0.2 for x = (•) 
X ~ 0.5 (o) and x = 1 (^)- Dashed lines represent the corresponding analytical expression 
from theorem\^and the solid line represents the function T{q) as defined in Eq. ((T]). is 
estimated from the average over 130 trials of 2-'^ cascades samples. Error bars are reported 
on the X = curve as vertical solid bar. These errors are of order of symbol size. 




Fig 3. Estimates of the function r^ip) of log-Poisson for x = (•) and x = 1 (°)- In 
both models we chose T = 2^'' and = 0.2. (a) Log-Poisson case with 5 = —0.1. (b) 
Log-Gamma case with (3 = 10. Solid lines represent the curves t{p). 
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X = curve. We can see that these errors are smaher or close to the symbol 
tickness. 

In Fig.Oare reported estimates of r^(p), x = 0, 1 for log-Poisson (fig.[3||a)) 
and log-Gamma (fig.[3|^b)) samples. The solide lines represent the theoretical 
t{p) functions for both models as provided by Eqs. (j45]) and (|46l) . We used 
the same estimation procedure as for the log-normal case. One sees that, 
in both cases, since the intermittency coefficient is the same for the three 
models, the classical To{p) curves are very similar to the log- normal curve 
(Fig. [2]). However, these models behave very differently in mixed regime: for 
X = 1, log-Poisson and log-Gamma both estimated scaling exponents be- 
come closer to the respective values of r(p) and are very easy to distinguish. 
Let us mention that such a analysis has been recently performed by two of 
us in order to distinguish two popular log-normal and log-Poisson models 
for spatial fluctuations of energy dissipation in fully developed turbulence 

m. 



APPENDIX A: LEMMA USED FOR THE PROOF OF THEOREM ?? 
Lemma 1. We have 



Mf\p) 



where L]^*" is defined by (125p and C is a constant that depends only on e. 



Proof. According to [25], if e S [0,1] and if {Xj}i<j<p are centered 
independent random variables one has 





p 


1+e- 




E 


i=l 




<cf]E[|X,|i+^' 

1=1 



where C is a constant that depends only on e (and neither on the law of 
X nor on P). Applying it with P = Nt = [2^^\ to the expression ()23p of 
Mj [p) , and using the fact that the random variables { [p] }m defined by 
([22]) are i.i.d, one gets 



< C2^^ E 



Mf\p)-E[M';"{p) 



Mf\p) 



,(0) 

j 

+ E 



Mf\p) 



Using the Jensen inequality we get the result. 



□ 
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Lemma 2. Assume that E < oo. Then we have for all m, 



E 



k=0 



where C is a constant that depends only on p and e. 

Proof. The proof of this result is very much inspired from [23]. Since the 
law of A^j™^ {p) is independent of m, we forget the supscript m throughout 
the proof. Using the definition (j22p . one gets 



(47) 

Let 
(48) 

then 



E 



\M,{p)\ 



l+e 



2-iP{i+^)E 



-•e{o,ip «=i 



r6{0,lp «=1 



rie{0,l}J r2e{0,l}J i=l 

It can be rewritten as 

(49) = Y + D, 

where Y corresponds to the non diagonal terms : 

riG{0,lp r2G{0,lp *=1 



(50) 



and D to the diagonal terms 
(51) 



rG{0,lp «=1 



The left hand side of (j47p is nothing but E 



. By writing that E 



\X 



E 



2\ 2 



using the sub-additivity of x i— > x^^^*^)/^, we get 



(52) 



E 



\X 



l+e 



< E 



, 1 + s 



+ E 



l + e 

D— 
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Let us first work with the Y term. We factorize the common beginning of 
the words ri and r2 in the expression (j50p of Y 

/ ^ / ^ 11 r\i / ^ 11 rri|i rr2|« 

A:=OrG{0,l}'= «=1 ri,r2G{0,lp"'= '=^+1 

ri|0^r2|0 

Again by the sub-additivity of x i— > x^^^*^)/^ and using the fact that the Wr\i 
are i.i.d., one gets 



E 



Y — 



fe=0 



(1+^) 



E 



/ ^ 11 rri|« rr2|«l 

'•l,'-2G{0,lp"* i=k+l 



{l+^)/2 



and by using Jensen inequality 



E 



fc=0 



(53) 



The variab^ 
the term E 



( E n ^ 

n,f2G{0,lp-'= «=A:+1 
■ri|07^r2|0 



^ E 

rGiO.ljfe 



E 



|^(n)^(r-2)| 



(l+^)/2 



es and are independant with finite expectation, thus 



|^(ri)2'(r2)| 



E we deduce: 

E 



is bounded by a constant C. Using 0^=^+1 ^ 



rr\\i rr2\i 



E[WP] 



(j-m+e) 



k=0 

E I E 

ri|0^r'2|0 

There are 2^ possible values for r and less than 2'^^^~^'^ values for the couple 
(ri,r2), thus 



E 



A:=0 
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Since 2-J^(p) = 2~^P2m [W^y 



(54) 



E 



\Y\~ 



k=0 



Let us now take care of the diagonal terms of X (jSip . First, we write 



rG{0,l}J «=1 



and using the E 



\Z 



l+e 



< oo we deduce that 



(55) E 

Merging ([Ml) and ([55]) into ([52]) leads to 



E 



\X 



fe=0 



and since E 



\MM\ 



E 



, it completes the proof. 



□ 



APPENDIX B: LEMMA USED FOR THE PROOF OF THEOREM ?? 

Lemma 3. Let / : [0, 1] — > M be some Borel function whose support is 

r_ (r+ 
2" ' 2° 



included in la^r = [^i -^^^^l for r € {0, 1}", a > 0. T/ien 



(56) 



(Moo,/> = 2-Mni^,|j(7lM,/ 



where f{x) = f{2 '^(x+r)) and //^^ zs a cascade measure on [0, 1] measurable 
with respect to the sigma field a{Wrr',r' G {0, 1}" ,a' > 1}. 

Proof. The scaling relation (j56p is easily obtained, by the definition of 
the measure Uoo , if f is the characteristic function of some interval / , , — j 

f*u(j) J a+a,rr' 

where r' G {0, 1}" , a' > 0. This relation extends to any Borel function / by 
standard arguments of measure theory. □ 

Lemma 4. Let /i : [0, 1] ^ M &e a piecewise continuous, non zero, func- 
tion. Then E [| (/Uoo, h) |^] > for all p > 0. 



imsart-aap ver. 2007/12/10 file: paper_ver5-3_IMS.tex date: May 2, 2008 



31 



Proof. By contradiction, assume that for somep > 0, E [j(/ioo) — 0- 
Hence (/^oo) h) = 0, P-almost surely. But using Lemma [3l 

= {floo, h) = (^floo, /il[0,l/2] ) + (^oo, /ll(l/2,l] ) 

where = h{2~'^-), = h{2-\- + l)) and/iS^/I^^ are independent 

cascade measures on [0, 1] . Thus we deduce Wq(^oo , 

almost surely, and since W > this shows that with probability one the two 
independent variables ('p'^,h^^^\ and (jl^ca , h^^'^\ vanish simultaneously. 



This is only possible either, if they both vanish on a set of full probabil- 
ity, or if they both vanish on a negligible set. Assume the latter, then the 
following identity holds almost surely 



Wi 



where the variables on right and left hand side are independent. These vari- 
ables must be constant, which is excluded by the assumption ¥{W = 1) < 1 
(recah (P). 

Thus we deduce that the variables (jl'^ , /i^*-* ^ are almost surely equal to 
zero. Hence: 



E 



0, for i = OA. 



Iterating the argument we deduce the following property: for any j > and 
k <2^ — 1, if we define a function on [0, 1] by h^^'''\x) = h{2~^{x + k)) we 
have 



E 



0. 



This is clearly impossible if we choose j, k such that h remains positive (or 
negative) on [k2~^ , (k + 1)2"-']. By the assumptions on h one can find such 
an interval, yielding to a contradiction. 

Lemma 5. We have 



□ 



E 



< C2^^E 



1+e 



where C is a constant that depends only on e. 

Proof. The proof is the exact same proof as for Lemma [TJ 



□ 
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Lemma 6. For any e > small enough we have, 

3 

A;=0 

where K is a constant that depends only on p and e. 

Proof. The proof basically follows the same lines as the proof of Lemma 
[21 The only difficulty, compared to this latter proof, comes from the fact that 
the quantity (/iooi5'j,fc) ^ priori involves several nodes of level j of the M- 
cascade. We have to reorganize the sum (jl3]). 

Since we are interested in the limit j — > +oo, we can suppose, with no 
loss of generality that j > J. In the following, we note 1^") the n-uplet 

^(n) = XI ... 1, where the 1 is repeated n times. 

The partition function (j43]) can be written Sf^^g{j,p) = J2k\(^f^oo-, djjjl^ i 
where the sum is over k € {0, !}•' such that ki = for some I < j — J- The 
sum can be regrouped in the following way, where a + 1 denotes the position 
of the last in the j — J first components of k: 

3-J-^ p 

S^,g{j,p)= E E E • 

a=0 rG{0,l}" sG{0,l}J 
q=rO 

We set 

q=rQ 

and consequently 

Actually, a exactly corresponds to the level of the "highest" node that is 
common for dyadic intervals in the support of g. Indeed, let us 

prove that 

(58) a > 0, Vs e {0, 1}^ Supp g.-^^^-j^^^ C Ia,r, 
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where q = rO. Indeed, according to (l42]) . the support of 5^. jrir^ is 
included in [2-JglO-"'"i-«)s, qlii-J-^-^) s + 2-(j--^)]. Then, 

i-J 

-- 2-'^r+ ^ 2-^ + 2"% 

i=a+2 

= 2^"r + 2^""^ - 2^(-''~'^) + 2"%. 



Since s varies in [0, 2"' — 1], and a < j — J — 1, it is easy to show that 



< 2~"-^ - 2-(J--^) 



and 



+ 2-^5, 
""^ + 2-% < 2"^ 



which proves (fSHI) . 

We are now ready to compute the upper bound for p)||j;^i+E(p 
Ea=o~"^ EsG{o,i}'^ \\Xs,a\\Li+^(¥)- Using ([57]), ([58]) and LemmaEl we get 



< 



re{0,l}" \i=l I 



where the [xQ, are independent cascade measures on [0,T]. 

Let us identify Xa^a with X as defined in Lemma [2] by ()48p in which j 



plays the role of a and Z^**) of 



. As in (HOl). we can 



^j-a,0lO--'-i-'')f 

decompose as the sum of the non diagonal terms Ya,s and the diagonal 
terms Da,s 

(59) =ya,s + ^a,s. 

Using the exact same development as the one we used for E [X^"'"'^] starting 
at Eq. (l47|l . we get the bound of the non diagonal terms corresponding to 



53)1 in which the term E 



_2-(rri)_^(rr2) 



has to be replaced by 



E 



(rri,s) 



9, 



j-a,Ol(-'-'^-i-"'; 



(rr2,s) 



9, 



which can be bounded (using (HH)) by A'2~^(-'~")('^(p)+^) . Going on with 
the same arguments as in Lemma [2l we finally get, the bound for the non 
diagonal terms corresponding to (|54|) 

a-l 

E 



I -'1,5 1 



i-1 

< ii'2"-^'(^+')^('')2(''"-'')(^+') 2-'=^(p(^+^))2''(^+')^(P^ 

fc=o 
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Following the arguments in Lemma [2] for the diagonal terms, we get 



(60) 



E 



Da 



1+1 

2 



By (I59])-(1601), we finally get 



E 



a,s I 



Then we write 



E 



*;=0 



a=0 se{0,l}-' 



and the lemma follows. 



l+e 



□ 



[1 

[2; 

[3: 

[4; 

[s: 

[6: 
[7; 

[9 

[lo; 

[11 
[12; 

[is; 

[14 
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